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We calculate the efi'ect of a quadratic term in the energy dispersion on the low-energy 
behavior of the Green's function of the spinless Tomonaga-Luttinger model (TLM). As- 
suming that for small wave- vectors q = k — kp the fermionic excitation energy relative 
to the Fermi energy is vpq + q'^ / {2m), we explicitly calculate the single-particle Green's 
function for finite but small values of A = qc/{2kp). Here kp is the Fermi wave- vector, 
(jc is the maximal momentum transfered by the interaction, and vp = kp/m is the 
Fermi velocity. Assuming equal forward scattering couplings 92 = 94i we find that the 
dominant effect of the quadratic term in the energy dispersion is a renormalization of 
the anomalous dimension. In particular, at weak coupling the anomalous dimension is 
7 = 7(1 — 2A^7), where 7 is the anomalous dimension of the TLM. We also show how to 
treat the change of the chemical potential due to the interactions within the functional 
bosonization approach in arbitrary dimensions. 

1. Introduction 

The exactly solvable Tomonaga-Luttinger model (TLM) has been extremely useful 



ill 



in the normal 
In this model 



to gain a better understanding of the electron-electron int, 
metallic state of electrons in one spatial dimension {d — 1 
only the forward scattering part of the electron-electron interaction is retained, 
which in d = Lean be parameterized by two phenomenological coupling constants 
(72 and (?4(RefMI). The second important feature of the TLM, which is crucial for 
its exact solubility via bosonizationErBu or other mcthodsQIa, is the fact that the 
energy dispersion of the electrons is exactly linear for all wave- vectors. Because the 
low-energy and long-wavelength physics is determined by the degrees of freedom in 
the vicinity of the Fermi points, it is reasonable to expect that the replacement of 
a general energy dispersion by the linear term in the expansion close to the Fermi 
points is justified. 

The single-particle Green's function GLx^ t) of the TLM in the space-time domain 
has been calculated by many authorsifltlDli. A detailed discussion of the rather 
complicated behavior of the Fourier transform G{k, lo) of G{x, t) in the momentum- 
frequency domain can be found in RefJ§. The structure of G{k,u}) is fundamentally 
different from the well-known behavior of the Green's function of a Fermi liquid: 
There is no quasi-particle peak, but G{k,uj) contains instead algebraic singularities, 
which can be interpreted physically in terms of bosonic charge- and spin- excitations 
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propagating with different velocities. The absence of a quasiparticle peak imphes 
also that the density of states of the TLM vanishes at the Fermi energy, and that 
the momentum distribution function exhibits only an algebraic singularity at the 
Fermi surface, in contrast to the jump discontinuity of a Fermi liquid. 

Haldaneli has emphasized that the above features are not special to the TLM, 
but are generic to the normal metallic state of interacting electrons in d = 1. He 
proposed to call such a state a Luttinger liquid, in contrast to the well-known 
Fermi liquid, which is the normal metallic state of interacting electrons in d = 3. 
In the seminal workd Haldane also emphasized that the various deviations from 
the idealized TLM can in principle be treated perturbatively within the framework 
of bosonization, and should only lead to a re-definition of the phenomenological 
parameters that appear in the Hamiltonian of the TLM, namely the Fermi velocity 
vp and the forward scattering couplings 52 and 54. For example, in bosonic language 
the non-linear terms in the expansion of the energy dispersion close to the Fermi 
points give rise to interactions between the Tomonaga-Luttinger bosons, which can 
be treated perturbatively within conventional bosonic many-body theoryS. In Ref.i 
Haldane performed an expansion of the equal-time Green's function to lowest order 
in the curvature parameter l/m, and pointed out that the non-linear terms in the 
energy dispersion can be expected to modify the anomalous dimension of tha TLM. 
However, the precise value of this renormalization was not obtained in RefJ3. This 
will be done in the pre sent work with the help of the functional integral formulation 
of bosonizationMll3ll3. 



2. Interacting electrons v^rith dominant forward scattering 



In this section we shall briefly describe how the single-particle Green's functio: 




G'(k, w) of interacting fermions with dominant Jopyaid scattering can be calculated 
within the functional bosonization formalismllil ijij in arbitrary dimensions. In 
particular, we shall rely on the results of Refs.lllliij, where a systematic method 
has been developed to handle the non-linear terms in the expansion of the energy 
dispersion close to the Fermi surface. It turns out for our calculation it is important 
to take the renormalization of the chemical potential due to the intera ctioii into 
account. This renormalization has not been discussed in previous worksy'E3E3li3'E3. 
Although we are ultimately interested in d = 1 , all equations in this section are valid 
in arbitrary d. 



2.1. Calculation of the single-particle Green's function via functional 
bosonization 

For simplicity, we shall ignore the spin degree of freedom. We also assume that 
the bare energy dispersion is Ck = and that the Fourier transform /q of the 
bare interaction between two electrons with momenta k and k' depends only on the 
absolute value of the momentum transfer k — k'. 

It is convenient to consider the Matsubara Green's function G(k, zu)„), where 
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ujn — 27r(n + ^)//3 are fermionic Matsubara frequencies. The inverse temperature 
is denoted by (3. Suppose that we are interested in G(k, at some momentum k 
in the vicinity of the Fermi surface. It is then convenient to measure all momenta 
relative to some reference point k" that is close to the Fermi surface, such that the 
length of q = k — k" is small compared with |k"|. Let us therefore define 

G"(q,ii;„)=G(k"+q,zc:>„), (1) 

where the index a indicates that q is measured relative to the reference point k" in 
the vicinity of the Fermi surface. The precise choice of k" will be discussed shortly. 
Note that the Fermi surface is defined as the surface in momentum space where the 
momentum distribution ^ 



n(k) - -_^G(k,zci„) (2) 



exhibits some kind of non-analytic behavior. Thus, in general we have to solve the 
many-body problem in order to know the precise location of the Fermi surface. Of 
course, if Luttinger's theorerrO'tS turns out to be valid in our system, then the 
volume of the Fermi surface is not changed when the interaction is switched on at 
constant density. For a spherically symmetric Fermi surface this implies that its 
radius is precisely identical with the Fermi wave- vector fc^ of the non-interacting 
system, 

kp = \/2m^o , (3) 

where /Xq is the chemical potential of the system without interactions at the same 
density. In general, however, the precise location of the Fermi surface is a priori 
not known. In fact, we cannot even exclude the possibility that for some type of 
interaction the momentum distribution is an analytic function of k, in which case 
the system does not have a Fermi surface. Note that in the shifted coordinate 
system the energy dispersion is 

ek"+q = ek-+v"-q-|-;^, (4) 
2m 

where v" = k"/m. Assuming that the system is confined to a finite volume V, the 
real-space imaginary-time Fourier transform of G"(q, icun) is 

G"(r, r)^^J2 e*'''-""""""^G"(q, ztD„) . (5) 

qn 

After the usual Hubbard-Stratonovich transformationS G*"(r,T) can be written as 
G"(r-r',r-r') = (a"(r,r',T,r')) , (6) 



where 



P2 

^" -Pr- -ek= -f /i-i0(r,T) 
Zm 



g"(r,r',T,T') = 5(r-r')r(r-r') . 

(7) 
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(a) (b) 

Fig. 1. Feynman diagrams which renormaUze the chemical potentiah Here the wavy-line is the 
screened interaction and the arrows represent the Green's function, (a) "Tadpole"-diagram; to first 
order in the interaction this is just the Hartree correction to the self-energy, (b) "Sunrise"-diagram; 
to first order this is the Fock correction. 



Here Pr = — iVr is the momentum operator, S*{t) — i g*'^"'^ is the antiperiodic 
(5-function, and (. . .) denotes functional averaging with respect to the effective action 
of the Hubbard-Stratonovich field (f>{r, r). The chemical potential of the interacting 
system is denoted by fj,. Let us now choose the reference point k" such that 

ek« -M-(i0(r,r)) . (8) 

Diagrammatically the term {i(j>{r, r)) corresponds to the sum of all tadpole diagrams 
which renormalize the chemical potential, see Fig.^a). The leading diagram of this 
type is just the usual Hartree diagram. We shall include the sum of these diagrams 
into the definition of the zeroth order Green's function. 



G'o(q, iw„ 



iujn - ek°-fq + M - («0(r, t)) 
1 



V" • q - 1^ 

1 2m 



(9) 



2m 

The effective action of the Hubbard-Stratonovich field is given by 



^off{0} = ^^/q-V-q,-»</'q,»-Trln[l-*Go(^] , (10) 

qn 

where the trace is over momentum-frequency space, with the infinite matrices Go 
and (f) given by 

[Go]qn,q'ii' — 5q,q' 5,i,n' Gq (q, itl;„) , (11) 
['^jqrijq'n' — 0q—q', n — n' ■ (-^^) 

Here 

Kn = / * / dre-^^-i— "-)[0(r,T) - (</.(r,T))] , (13) 

py J J-13/2 

where a;„ = 27rn//3 are bosonic Matsubara frequencies. In Eq.([l3|) we have sub- 
tracted the average ((/)(r, r)), which by translational invariance is independent of r 
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and T. This subtraction modifies only the zeroth Fourier component 0o = 0q=o,ri=o 
such that {(f>o) ~ 0. Note that for non-zero q or ujn translational invariance implies 
= 0. 

FromEqs.(|l|) and (|) it is clear that the Fourier transform G(r, r) of the physical 
Matsubara Green's function G(k, ia)„) is simply 

G(r,r) = e''^°-"-G"(r,T) . (14) 

If one works with locally linearized energy dispersion the Fermi surface has to be 
subdivided into several sufficiently small sectorsllirEj, which are labelled by the 
discrete index a. In this case the right-hand side of Eq.(14) should be summed 
over all these sectors in order to recover the physical Green's function. Note that 
the sectorization would require the introduction of cutoffs in momentum space. 
However, for our purpose such a subdivision into sectors is really not necessary as 
long as we restrict ourselves to the calculation o/ G(k° -t-q, w) for |q| <C |k"|. Of 
course, without sectorization the curvature of the Fermi surface has to be taken into 
accounty-ti. 

It turns out that for interactions that involve only small momentum transfers 
|q| <C |k"| it i s s iifficiput to perform the averaging in Eq.(H) within the Gaussian 
approximationliinia, which amounts to an expansion of the logarithm in Eq.(|lO|) to 
the second order. Within this approximation the effective action is given by 

-Seff ~ iPVpo4>0 + ^ X! ■^RPA(q' iw„)0_q,_n</)q,„ , (15) 

qn 

where 

Po = ^Y.^oin,i^n), (16) 

qn 

/RPA(q,iw„) = — — A — . — ^ (17) 
1 + /qno(q,?t^„) 

is the effective interaction within the random- phase approximation (RPA). The 
"non-interacting" polarization can be written as 

no(q, iLU„) = E ^0 (q'. i^n')G^{^' + q, i^n'+n) ■ (18) 

q'n' 

Let us now focus on the solution of Eq.(0). Making the generalized Schwinger- 



and 



g"(r, r', r, r') = G^r, r', r, r')e^"<^r,r)-^-ir' .r') ^ (^g) 
we obtain with k" given by Eq.(p|) 



Zm m 



er(r,r',r,T')+e?(r,r',r,r') 
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-dr - V" • Pr - 

5{r-v')5*{T -t') 



p2 
r 

2m 



$"(r,T)-z[0(r,r)-(0(r,T))] 



[A"(r,T) 
2m 



Here 



A"(r,T) = Pr<i>"(r,r) 



(20) 



(21) 



is a longitudinal vector potential. Eq.(19) is a gauge transformation. Obviously, 
there are infinitely many different choices for the gauge factor $"(r,T). For exam- 
ple, if the bare interaction /q is the three-dimensional Coulomb interaction, then 
A" (r, r) can be identified with the longitudinal part of the usual ^ctor potential 
of electromagnetism if we choose — 9T-$"(r, r) — i0(r, r). In Ref.c3 it was shown 
that this choice is useful to resum the leading logarithmic singularities in the per- 
turbative expansion of the Green's function of two-dimensional disordered electrons 
that interact with long-range Coulomb forces. To include the renormalization of the 
chemical potential into the functional bosonization approach, a slight modification 
of the gauge choice given in Refs.EUll is necessary. It is easy to see that Eq.(po|) 
can be satisfied if 



-dr 



2m m 



- Do 



5{r-r')5*{T-T') , 



(22) 



Pi 

2m 



$"(r, r) 



= i[0(r,r)-(0(r,T))-0o] + 

Here (f)o is the zeroth Fourier component of 0(r, 
the average diamagnetic energy. 



[Pr$"(r,r)]2 



Do 



2m 

((/>(r, r)) (see Eq. 



Do = 



1 

w 

I 

pv 



dr 



= dv 



[P,ci>"(r,r)]^ 
dr 

-p/2 2m 

^ JA"(r,r)]^ 

-p/2 2m 



(23) 



and Do 



(24) 



By construction the Fourier component corresponding to q = and w„ = of the 
right-hand side of Eq.(p3[) vanishes, so that in the iterative solution of the eikonal 
equation ( p3[ ) ambiguities due to vanishing eigenvalues of the differential operator 
on the left-hand side are avoided. Such an e xplicit separation of the zeroth Fourier 
component, which was not discussed in RefsO'E3, is important to incorporate the 
renormalization of the chemical potential into functional bosonization. 

At this point we solve the eikonal equation ( ^^ to first order in (j){r,T), which 
is sufficient to calculate the Debye- Waller factor arising after the averaging from 
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the gauge-factor in Eq.(p^) to first order in the RPA interaction. The functional 
averaging of the "prefactor" Green's function is then perfonned within the self- 
consistent Born approximation, as described in detail in Ref.Ell. The final result for 
the G reen 's function of the interacting many-body system for wave- vectors close to 
k" is 00 



G(k°+q,icD„) = dr dre~'('>-"-'^"")G"(r,T)e'3? , (25) 

J J-B/2 



'-/3/2 

with the Debye- Waller factor given by 



s 1 /R,PA(q,^cJ„)[l-cos(q■r-a;„T)] 

/5^tr[-""V".q-|i][z..„-v".q+^] ' 

and the prefactor Green's function 



qn 



Here 



(q-q')q'^ ^ ( <i q.' \2 



(26) 



G"(r,T) = ^^e^('i-"--"^)[l-Kr"(q,*ci„)]G?(q,z(:-„) . (27) 



Gi (q, ^^n) = , , (28) 

S^(q,iw„) = -^^yRPA(q',iw„')G'i(q + q',iw„+„') 



(29) 



F"(q,i(D„) = ■^^/RPA(q',iw„,)G^(q-f q',i(D„+„0 

q'n' 



3--f 2 



qq 



and 



q'-a2][iw„--v"-q' + a:i] 



= (-Do) = -^X!'^RPA(q,«w„) 



qn 



q 

2m 



(30) 



(31) 



Note that 5^ p can be identified with the renormalization of the chemical potential 
arising from the usual "sunrise" -diagram shown in Fig.|l|(b). To leading order in the 
bare interaction, this is just the Fock diagram. Note that the Hartree renormaliza- 
tion of the chemical potential has already been taken into account in the definition 



8 Tom Busche and Peter Kopietz 



of the zcroth order Green's function by a suitable choice of the reference point k", 
see Eq.(||). This different treatment of the two renormahzations of the chemical 
potential is justified if we formally treat 1/m as a small parameter. Obviously, 
Slip oc 1/m, whereas for the Hartree correction that appears in Eq.(^) we obtain to 
leading order 

S^IH ^ (th) ^ PofRPA{0,0) ■ (32) 

Keeping in mind that in d dimensions po oc kp (x m'^, it is clear that the Hartree 
renormalization of the chemical potential diverges if we let m ^ oo holding vp = 
kp/m constant. This is the reason why this renormalization cannot be treated 
within a l/m-expansion, but should be included in the definition of the zeroth 
order Green's functioncll. 

It is important to emphasize that the above exnansion is quite different from the 
expansion in powers of \/m proposed by Haldanea. The advantage of our method 
is that in the absence of interactions we recover the exact free Green's function 
with non-linear energy dispersion, which contains of course infinite orders in 1/m. 
In Sec. 3 it will become clear that the effective small parameter in our expansion is 
proportional to /rpa/^i, so that we obtain the exact result in both limits /rpa 
or 1/m ^ 0. 

2.2. Luttinger's theorem 

The Fermi surface of an interacting Fermi system can be defined as the surface in 
momentum space where the momentum distribution defined in Eq. (^ has some kind 
of non-analyticity in the zero-temperature limit /3 — > oo. The Fermi surface is in 
general not identical with the surface defined by Eq.(^), so that our reference point 
k" is not located precisely on the Fermi surface. Note, however, that according 
to Luttinger's theoremlill the volume enclosed by the Fermi surface is not changed 
by the interaction, so that the radius of a spherically symmetric Fermi surface is 
simply kp — \/2m/io- On the other hand, our reference point k" is defined by 
Ek" = M ^ 5p,H^ so that to leading order 

|k"| - kp ^ 6hh - fJ.o)/vp = Sfip/vp , (33) 

were we have used the fact that to first order in the RPA screened interaction 
= /io + 6fiH + 6fip, see Eqs.(^l]j3|). It is now easy so see that for sufficiently 
regular interactions the singularity in the momentum distribution occurs at the 
non-interacting kp. Substituting Eq.(|25|) into the definition (^ we obtain 

n(k" + q) = y" dre-*i ""G"(r, 0)e'^"^''-°'^ . (34) 

Obviously the behavior of this function for small q is determined by the long- 
distance behavior of the prefactor Green's function G"(r, 0) and of the Debye- 
Waller factor Qiix, 0). In a Fermi liquid (5i(r, 0) approaches a finite constant at 
large distanceal^'E2l, while for a Luttinger liquid the Debye- Waller factor exhibits 
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logarithmic singularities, which wiU be discussed in Sec. 3. 3. The long-distance 
behavior of G"(r, 0) is dominated by the singularities of its Fourier transform given 
in Eqs.(p7[) and (p^). The important point is now that for finite S/ip the singularity 
occurs at a finite value of q. This is of course due to the fact that we have not 
measured the momenta relative to the true Fermi surface. It is therefore convenient 
to shift the momentum q in Eq.(p8[) by setting 

q = q' - (5k" , (35) 

where i5k" should be chosen such that it cancels the chemical potential shift in 
Eq.(p8|). For our spherical Fermi surface this means that 

^^|k"| - VF |k"| • ^-^^^ 

Renaming again q' — q, we obtain 

G"(r,0) = e-*''''" "-G^(r,0) , (37) 

with 



qn 

l + r"(q,iw„) 



»c:;„-(v"-^)-q-^-S?(q,z(:>„ 



(38) 



Since we retain only terms to first order in /rpa, we have ignored the shift 5k" in 
the arguments of and F". This is justified, because according to Eqs.( |3l|j3^ ) 
the shift (5k" is already of order /rpa- The momentum distribution can then be 
written as 

n(k" - (5k" + q) = y" (ire~*i ""G^(r, 0)e^?('''°) . (39) 

The important point is now that for a conventional Fermi liquid the right-hand 
side of Eq.(^9|) is non- analytic at q = 0. Because by construction |k" — (5k" | = 
kp — y/2/rnJIo, the momentum distribution is non-analytic at the non-interacting 
Fermi surface, in agreement with Luttinger's theoremHI. In Sec. 3 we shall explicitly 
evaluate the above expression in one dimension, and verify that even in this case 
Luttinger's theorem remains validtSl. 

The reader may find our discussion in this section rather pedantic. However, as 
will become evident in Sec. 3, for the calculation of the effect of the quadratic term 
in the energy dispersion on the spectral function of the TLM it is crucial to take the 
shift of the chemical potential due to the interaction into account. Of course, if the 
location of the true Fermi surface is known a priori (this is the case, for example, 
in systems where Luttinger's theorem is valid and the Fermi surface is known to 
be spherically symmetric), it seems more natural to choose a point on the true 
Fermi surface as the reference point for the linearization of the energy dispersion. 
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However, for general interactions and energy dispersions the location of the Fermi 
surface is not known a priori. In this case the general method outlined above can 
be used to calculate the true Fermi surface perturbatively. In the following section 
we shall explicitly demonstrate how this method works in one dimension. 



3. Non-linear energy dispersion in one dimension 

In d = 1 there are only two Fermi points, which can be labelled by a = ±. The 
zeroth order Green's function defined in Eq.(H) can then be written as 

(q, ^^n) = -T- ■ (40) 

lujn ~ avpq - ^ 

The reference points fc" and the renormalized Fermi velocity vp are (see Eq.(p|)) 



A;" = Q!\/2m(/io + (J^f) , 

VF 



(41) 
(42) 



^2(^0 + 5^F)/m . 

It is important to emphasize that according to the generalized Luttingcr theoremEl 
even in d = 1 the volume of the Fermi surface is invariant as we switch on the 
interaction at constant density. Hence, the Fermi wave- vector kp is still given 
in terms of the chemical potential of the non- interacting system, kp = \/2mfiQ. 
Therefore the velocity vp — \k°'\/m is not identical with the bare Fermi velocity 
Vp = kp/m. 

In this section we shall explicitly evaluate the expressions of Sec. 2 in one dimen- 
sion for a bare interaction of the form 



fq /o0(9c - \q\) , qc<t:kp 



(43) 



Note that this corresponds to the TLM with forward scattering interactions g2 {q) = 
gi{q) = fq- In one dimension the polarization function in Eq.([l^) can for V ^ oo 
and /? ^ cxD be written as 



Ho((j,icj) = - — In 
2'Kq 

In the limit 1/m ^ this reduces to 

1 



[vpq 



q \ 

2ml 



{vpq ^ 
vpqf 



q \ 

2m > 



TTVp (vpq)^ + UJ^ 



1/m . 



For finite m the RPA interaction can be written as 

TTVpFq 



/rpa('7, j^) = 



where we have introduced the dimensionless interaction 

A 

TTVp 



Fg^Foeiq,-\q\) , Fo = ^ 



(44) 



(45) 



(46) 



(47) 
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and the dimensionless polarization 



P(k, iu) = — In 
4fc 



(48) 



3.1. Chemical potential 

Let us first consider the Fock-renormahzation of the chemical potential given in 
Eq.(^l]), which turns out to be finite even in d = 1. Introducing the dimensionless 
variables p = q/ qc, u = lj / {vpq) and the small parameter 



\ = 7r^, (49) 

Zmvp 



we obtain from Eq.(|T]) 



where 



duF = -Do) = / dp\p\ / du- ^— ^ — — , 50) 

87rm7_i 1 + FoP{Xp,iu) 



A{k, lu) = — — . (51) 

[lu — a — k\ [lu — a -\- k\ 



For small A we may expand the functions P(Xp,iu) and A(Xp,iu) in powers of A. 
To leading order we obtain 

P{\p,iu) « (52) 

A(Ap,m) « ^ (53) 

[tu — a) 

Recall that a — ±1 labels the two Fermi points. The integrations in Eq.(^0|) are 
then easily performed, and we obtain 

where 

7 = ^ (55) 

2v/l + i^o(Vl + -Fo + l)2 

will be identified below with the anomalous dimension of our model. Using po = 
mvp/'K the Hartree renormalization of the chemical potential can be written as 

SfiH {i4>o) = mw|- ^ , (56) 

1 + -Fo 

so that the total shift of the chemical potential is 
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Note that the Hartree renormahzation is divergent in the Tomonaga-Luttinger hmit 
TO — > cxD with kp/m — const. This is the reason why it is necessary to include this 
renormahzation into the definition of the zeroth order Green's function. Note also 
that by construction the prefactor Green's function ( p8| ) depends on the bare Fermi 
velocity vp — kp/m. This follows directly from Eqs.( p6| , p2| , p4[ ), which imply in 
d — \ for the renormalizcd Fermi velocity 

ip^vp[l + X^^], (58) 

and hence for the relevant velocity that appears in the prefactor Green's function 

a. . Sup 

V = avp — a = avp . (59) 

TO mvp 



Here 



A = (60) 
2mvp ^ ' 



is defined in terms of the bare Fermi velocity. Obviously 

A = A + 0(A3), (61) 
so that to leading order in an expansion in powers of A we may replace A — > A. 

3.2. Prefactor Green's function 

Next, we calculate the functions S]"(g, iw) and F"((7,iw) defined in Eqs.(p9|) and 
( pO| ) foT V oo and P oo. For small values of q and lu and to leading order in 
/rpa it is sufficient to replace G" — > Gq on the right-hand sides of Eqs.(p9|) and 
(pO|). This neglect of self-consistency is justified a posteriori: because the result of 
the lowest order Born-approximation is finite and small, we presume that the self- 
consistent Born-approximation is not necessary. Introducing the same dimensionless 
variables as above, we obtain 

2 



2-,i[q,iULi) — —vpq — / ap\p\ / du 

71- J-1 J-oo 



l + FoP{Xp,iu) 



P+f- 

He 



G^^ +p[iu-a-\{2f^+p)] 



(62) 



^ [ , r , FoA{Xp,iu) 



A f f 

Y°'{q,ioj) = - dp du 

TTJ-l J-oo 



1 + FoP(Ap, 



H+2sgn(p)^ 

(63) 



G^'+p[iu-a-Xi2f^+p)] 

where ^ 

. , iuj — avpq — 77— 



Q-l ^ 12 2m /g4) 

VfQc 
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To make further progress, we expand Eqs.(|62[) and ( |63| ) in powers of A. To leading 
order, we may simply set A = in the integrands and replace vp ~^ vp, A — *■ A. The 
integrations can then be performed analytically, but the result is rather complicated 
and not very illuminatingr3. For our purpose, we only need the leading terms for 
small q and lu, which are 



5vi = —avpX^'f- 



(65) 
(66) 



where 7 is the anomalous dimension of the TLM, which is obtained from Eq.(|55|) 
by replacing Fq ^ Fq = fa/iTTvp)- For the function Y{^{q,iuj) we obtain to leading 
order 

' 9c 



r"(g,ic.) = A 



ci- 



vpqc 



q 

C2 — 

9c 



(67) 



where ci and C2 are numerical constants that depend on Fq 



3.3. Debye- Waller factor 

Finally, consider the Debye- Waller factor defined in Eq. ( p6| ) . Introducing again the 
dimensionless integration variables p — q/qc and u = Lu/{vpq), we obtain in d = 1 

47ri_i|p|i_^ l + FoP{Xp,tu) 



X [1 — cos(p(i — uf))] 



(68) 



where x = qcX and f = vpqcT. Since we are only interested in the leading behavior 
of Q"{x,t) for large x or r, we may set A = in the integrand of Eq.(|68|). This is 
easily seen from the fact that the factor of \/\p\ in the integrand is responsible for 
a logarithmic growth of (5"(x,t) for large x or r. If we now expand the integrand 



of Eq.(68) in powers of A, we obtain additional powers of p, so that the resulting 
contributions to Q"(a;,r) are bounded and hence can be neglected as far as the 
leading logarithmic behavior is concerned. In the limit |r| ^ (t'F9c)^^ and |a;| 3> 



the leading term is 



Q?(a;,r)^^ln 



x"^ + {VcTY 



In 



X + iavpT 



x + laVcT 



where 



vpyl 



(69) 



(70) 



is the velocity of collective charge oscillations and 7 is the renormalized anomalous 
dimension, see Eq.(|55|). Eq.(|69|) the well-known result of the TLM, but with renor- 
malized parameters 7, Vc and vp. A more rigorous way to obtain Eq.(|69|) is to first 
perform the u-integration in Eq. (^) without expanding the integrand in powers of 
A. Then, in the same limit as above, namely for large x and t, we can also perform 
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the p-integration analytically, which leads precisely to the presented resultEl. Note 
that Eq.(69) implicitly depends on our small parameter A, because the renormalized 
parameters 7, Vc, and vp are A-dependent. Hence, Eq.(|69|) goes beyond a simple 
expansion in powers of A. 



3.4. Spectral function 

Combining all the results of the previous subsections, we finally obtain for the 
Matsubara Green's function for wave- vectors q close to the two Fermi points ±kp — 



dx 



X + lavpT 
X + iavcT 



dTe-'^'>''-'^''^G'^(x,T) 



with 



Here 



G^{x,t) 



dq 
2^ 



duj 
— ( 
27r 



X^ + (VcT)^ _ 

1 + Y"{q, iuj) 



c^v,q[^ + \^] 



Vi = Vp 



1 - A27- 



1 



O(A^) 



(71) 



(72) 



(73) 



where we have used Eqs.( p9| , |66| ). Because the anomalous dimension in Eq.(|69|) 
depends on the effective dimensionless coupling — fa/{T:vF), which involves the 
renormalized Fermi velocity vp , we obtain a small correction to the anomalous 
dimension 7 of the TLM. Keeping in mind that our result for the renormalization 
of the chemical potential is based on a weak coupling expansion, we obtain from 
Eq.(|69|) in the Tomonaga-Luttinger limit A ^ and for small Fq, 



P. 



Hence, for finite A the effective anomalous dimension is 



7 



7 



1 + 2A27 



7[1-2A27 + 0(A3)] 



(74) 



(75) 



This is the leading correction to the anomalous dimension of the TLM due to the 
quadratic term of the energy dispersion at weak_coupling. We expect that the non- 
Gaussian corrections to the Debye- Waller factorllil do not change the correction term 
in Eq. ([75|) : simple power counting shows that these non-Gaussian corrections do 
not contain any infrared divergences, so that they cannot give rise to logarithmically 
growing corrections to the Debye- Waller factor, which are necessary to renormalize 
the anomalous dimension. 

Because our result for the vertex function Y^^q^iLo) in Eq.(|67|) vanishes in the 
limit (7, w — > 0, it does not modify the qualitative behavior of the spectral function for 
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small A (Ref. B). To leading order, we may therefore set Y"{q, ioj) 0. Restricting 
ourselves to the limit of large x and r and neglecting the term of order in the 
denominator of Eq.(|7^), we obtain 

GU^,r) = ^—^ , \xl\r\^^. (76) 

The imaginary-time Green's function in Eq.([7l|) becomes then 

-1 /"OO />oo 

Giakp+q^ioj) = / dx dTe-*(«^-"^) 

ZTTJ 



1 



X + iaVcT 



oo 
-2 



x'^ + {VctY 



7/2 ^ - 

a; + lavpT 
X + iaviT 



(77) 



Note that for finite A and Fq all three velocities vi, vp and Vc are different. Only 
for linearized energy dispersion (i.e. for the TLM) vp — Vi, so that the last factor 
in Eq.([77|) is exactly unity. For finite A, however, the velocity degeneracy is broken, 
which has rather spectacular consequences for the spectral function: there appears 
an additional algebraic singularity in the spectral function, which is not present in 
the TLME3. The fact that even a small breaking of the velocity degeneracy in a 
special model of interacting fermions in d = 1 can lead to new singularities in the 
spectral function has recently been pointed out by Ho and Colemanil. Because 
in our model the curvature of the energy band implies a velocity dispersion, it is 
not unreasonable to expect that the curvature leads to new features in the spectral 
function, which are not present in the TLM. However, at this point we cannot 
exclude the possibility that this new feature, whose appearance depends in a rather 
subtle way on a small renormalization of effective velocities, is completely washed 
out by non-Gaussian corrections or by the subleading corrections to the Debye- 
Waller factor, which have been neglected in our leading orderxalculation. Moreover, 
this singularity has a very small weight (which can be shownEa to vanishes as A^ for 
A — > 0), so that it is not of any practical importance. 

For these reasons we shall replace the last factor of Eq.(|7^) by unity in the 
following analysis. Note, however, that the charge velocity Vc and the anomalous 
dimension 7 are defined in terms of the renormalized velocity vp — vp{l + A^7), 
which is slightly larger than the bare Fermi velocity vp — kp/m. To calculate the 
spectral function 

A°'{q,uj) =--lmG{akp+q,uj + iO+) (78) 

TT 

from Eq .(j77| ) with wj;^ = ui, it is convenient to first consider the real time Green's 
functioiJHJ 

G%{x,t) ^ie{t)G"{x,iT ^ -t + iO) . (79) 
The spectral function for w > is then given by 

G'^{q,Lu)^27re{u;)A''{q,Lu) , (80) 
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where 



We obtair 



dt 



(81) 



A'^{q,to) 



7 



2vcqcJ 2r2(l + 2) 
X - \vcq\){co + avcq)^ {uJ - avcq)^^^ 



(82) 



Note that 7 can be expanded in powers of A, but a perturbative calculation of 
A°'{q,u!) in powers of A would lead to unphysical logarithmic singularities. To 
illustrate this point, consider the momentum distribution. From Eq.(^4[) it is easy 
to show that 



n(akp + q) 



1 



Casgn((7) 



(83) 



where C is a numerical constant. Using 7 = 7 — 2\^^'^ , an expansion of the second 
term in powers of A yields 



(84) 



Thus, a naive expansion of n^akp + q) in powers of A would generate logarithmic 
terms, which become arbitrary large for sufficiently small q. Our method effectively 
resums all these corrections. 



_9 


1 


(]_ 


7 


l + 2A27ln 
















q 



4. Conclusion 

In this work we have studied how in one dimension the quadratic term in the expan- 
sion of the energy dispersion close to the Fermi points modifies the spectral function 
of the TLM. The most important effect of the non-linearity is a renormalization of 
the anomalous dimension, which we have explicitly calculated to leading order in 
the small dimensionless parameter A — qc/{2mvF)- Note that this implies that cor- 
relation functions can in general not be expanded in powers of the A, because they 
exhibit algebraic singularities. This leads to a non-analytic A-dependence. We have 
also pointed out that the non-linear terms in the energy dispersion might give rise to 
new features in the spectral function, which are not present in the spectral function 
of the TLM. However, we cannot exclude the possibility that these non-universal 
features are washed out by higher order corrections which have been neglected in 
our calculation. 

For a proper treatment of the quadratic term in the energy dispersion it is cru- 
cial to take the renormalization of the chemical potential due to the interaction into 
account. We have shown how to include this effect into the functional bosoniza- 
tion approachEJ in arbitrary dimension, emphasizing that the "tadpole" diagrams 
describing a Hartree-renormalization of the chemical potential have to be treated 
exactly. 
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